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ABSTRACT 

We run numerical simulations of molecular clouds (MCs), adopting properties similar 
to those found in the Central Molecular Zone (CMZ) of the Milky Way. For this, we 
employ the moving mesh code Arepo and perform simulations which account for a 
simplified treatment of time-dependent chemistry and the non-isothermal nature of 
gas and dust. We perform simulations using an initial density of no = 10^ cm“^ and a 
mass of 1.3 x 10^ Mq. Furthermore, we vary the virial parameter, defined as the ratio 
of kinetic and potential energy, a = ^kin/l^potl, by adjusting the velocity dispersion. 
We set it to G = 0.5, 2.0 and 8.0, in order to analyze the impact of the kinetic energy on 
our results. We account for the extreme conditions in the CMZ and increase both the 
interstellar radiation field (ISRF) and the cosmic-ray fiux (CRF) by a factor of 1000 
compared to the values found in the solar neighbourhood. We use the radiative transfer 
code Radmc-3d to compute synthetic images in various diagnostic lines. These are 
[Cll] at 158/im, [Ol] (145/im), [Ol] (63/im), ^^CO (J = 1 -> 0) and ^^CO (J = 1 ^ 0) 
at 2600 /rm and 2720 /im, respectively. When a is large, the turbulence disperses much 
of the gas in the cloud, reducing its mean density and allowing the ISRF to penetrate 
more deeply into the cloud’s interior. This significantly alters the chemical composition 
of the cloud, leading to the dissociation of a significant amount of the molecular gas. 
On the other hand, when a is small, the cloud remains compact, allowing more of 
the molecular gas to survive. We show that in each case the atomic tracers accurately 
reflect most of the physical properties of both the H 2 and the total gas of the cloud 
and that they provide a useful alternative to molecular lines when studying the ISM 
in the CMZ. 
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1 INTRODUCTION 


The inner few hundred parsecs of our Milky Way, known 
as the Central Molecular Zone (CMZ), are rich in molecular 
and dense gas and account for about 5% of the molecular gas 
content of the Galaxy ( Morris Serab}^|1996 ). As shown 
by Longmore et al. (2013), the CMZ contains gas with den¬ 
sities of the order of several 10^ cm“^ and a total mass of 


10^ Mq. Furthermore, it is highly turbulent (Shetty et al. 


2012) and illuminated by a strong interstellar radiation field 


(ISRF) combined with a high cosmic-ray flux (CRF) (Yusef- 
[Zadeh et al.|2007[ Clark et al.|2013 ). It thus provides an ex¬ 
cellent laboratory for studying the physics of the interstellar 
medium (ISM) under extreme conditions (such as those that 


also occur in distant starburst galaxies) with very high res¬ 
olution. 

Previous observations suggested the existence of copi- 
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velocity structure of the underlying molecular component of 
the ISM. However, the effect of the extreme physical condi¬ 
tions in the CMZ on the CO distribution within the gas and 
the degree to which CO is a biased tracer of the underlying 
cloud properties remain relatively unexplored issues. Pre¬ 


vious numerical models of clouds in the CMZ (Clark et al. 
2013) or in other harsh environments ( Clark Glover|20I5 ) 
have shown that the gas is likely to be highly chemically in- 
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homogeneous. This means that molecular tracers might tell 
us only little about the kinematics, the temperature and 
density distributions of the cloud and so on. Hence, in order 
to learn more about the internal physics (e.g. heating, cool¬ 
ing and the chemistry) of the warm gas, which makes up a 
large fraction of the total mass, we have to use other probes. 

Recently, [Clark et al. (2013) have shown that the [Ol] 
63 jim. and the [Cii] 158 /xm fine structure lines dominate the 
cooling in CMZ clouds over a wide range in densities. More¬ 
over 


which we use for computing synthetic observational maps of 
the clouds in various diagnostic lines. 

2.1 Hydrodynamical and chemical model 

We use the moving mesh code Arepo ( Springel|2010 ) in or¬ 
der to run numerical simulations of various model clouds. 
We have added to Arepo a detailed atomic and molecular 


cooling function, described in detail in Glover et al. (2010) 


of these fine structure lines in the CMZ, carried out with the 
ISO satellite. They show that both lines are very strong in 
this region and can be used to infer important information 
about the physics of the gas in the CMZ. Hence, using such 
atomic tracers in order to study the physics of MCs in the 
CMZ seems to be a promising approach. 

In this paper we follow up on this idea and investigate 
the applicability of atomic tracers to study the properties 
of the ISM in the CMZ. We aim to provide synthetic obser¬ 
vations, readily comparable to observations of fine structure 
lines. To do so, we compute synthetic observations of MCs 
in a CMZ-like environment using the radiative transfer code 
Radmc-3d ( |Dullemond| 2012| for various diagnostic lines. 


Rodrfguez-Fernandez et al. ( 2004 ) present observations and Glover & Clark (2012), and a simplified treatment of the 


specifically for [Cll] (158/xm), [Ol] (145/xm), [Ol] (63/xm), 
^^CO (2600/xm) and ^^CO (2720/xm). We perform simula¬ 
tions of various clouds with the moving mesh code Arepo 
( Springel||2010| ) using environmental properties similar to 


those experienced by a typical CMZ cloud. In particular, we 
adopt values for the ISRF strength and the CRF compara¬ 
ble to those inferred for the dense CMZ cloud GO.253+0.016, 


also known as “The Brick” (see e.g.|Gusten et al.||1981| |Lis 

et al. 

1994||Lis & Menten|19981|Lis et al.|2001||Molinari et al. 

2011 
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[Longmore et al.||2012| [Kauffmann 

et al.||2013| Clark et al. 

|2013| [Johnston et al.[|2014| Rath- 

borne et al.|2014| 2015). 

We use values for the ISRF and the 


CRF that mimic the harsh conditions assumed to be found 
in the CMZ, which are a factor of ~ 1000 larger than the 
values measured in the solar neighbourhood ( Habing||l968 


Draine|1978 Mathis et al.|198'3 ). We simulated clouds with 
an initial number density of no = 10^ cm“^ and studied the 
impact of different virial a parameters, a = 0.5, 2.0 and 8.0, 
on the physical properties. More information about the sim¬ 


ulations can be found in Bertram et al. (2015). Furthermore, 


we note that this study is the first of two papers which an¬ 
alyze the properties of CMZ-like clouds. In this paper, we 
focus on maps of column density and integrated intensity, 
while in Paper H (Bertram et ah, in prep.), we study the 
kinematic properties of our model clouds. 

This paper is structured as follows. In section we 
present our numerical simulations and the radiative trans¬ 
fer post-processing tool. In section we show and discuss 
the results of our studies. We present a summary and our 
conclusions in section |4] 


2 METHODS 

The numerical simulations studied in this paper are de¬ 


scribed in more detail in Bertram et al. (2015). However, we 


summarize the most important aspects of the simulations 
here and introduce the radiative transfer post-processing. 


basic chemistry of the gas. The chemical network is based on 


the work of Nelson & Langer (1997) and Glover & Mac Low 


(2007), and allows us to follow the formation and destruc¬ 


tion of H 2 and CO self-consistently within our simulations. 
The network tracks the abundances of 6 species and fol¬ 
lows 14 chemical reactions. The simplified [Nelson Langer| 
(1997) network is known to somewhat overestimate the rate 


at which CO forms, owing to its neglect of atomic carbon 
(see the detailed discussion in Clover k, Clark|2012 ). It also 
neglects the reaction 

CO + He+ C+ + O + He, (1) 

which plays an important role in regulating the CO abun¬ 


dance in gas exposed to a high CRF (see, e.g. Clark & Glover 
2015 Bisbas et al.||2015 ). Our simulations thereby overesti¬ 


mate the CO abundance in the simulated clouds, and hence 
potentially underestimate the C^ and O abundance in re¬ 
gions where the CO abundance is large. Therefore, in order 
to establish the effect that this has on the synthetic maps, 
we make two versions of each map: one in which we use val¬ 
ues for the fractional abundances of C^ and O taken from 
the simulations and a second version in which we assume 
that all of the carbon is in the form of C^ and all of the 
oxygen is present as O. This can be achieved by assum¬ 
ing zero CO abundances in our simulations. As discussed in 
more detail in Appendix [C] we hnd only minor differences 
in the [Oi] (145 /xm) emission PDF and negligible differences 
in the [Oi] (63/xm) and [Cli] (158/xm) PDFs. We therefore 
conclude that the known weaknesses of our simplified chem¬ 
ical model have little influence on the results we obtain for 
[Cll] and [Ol]. Full details of the chemical model with a de¬ 
scription of how the chemistry interacts with the ISRF via 


the TreeCol algorithm can be found in Clark et al. (2012) 


Further examples of the use of our chemical model with the 
Arepo code can be found in Smith et al.| (2014a) and Smith 
eT^ (2014b). 


We assume that the gas has a uniform solar metallicity 
and adopt the standard ratio of helium to hydrogen. The 
abundances of carbon and oxygen are taken from [Sembach] 


et al. 


(2000). We use = 1-4 x 10 and xo — 3.2 x 10 


where xc and xo are the fractional abundances by number 
of carbon and oxygen relative to hydrogen. We note that 
the CMZ has super-solar metallicity. Nevertheless, we use a 
uniform solar value in order to be conservative regarding the 
cooling rates in our runs. When we start the simulations, hy¬ 
drogen, helium and oxygen are in atomic form, while carbon 
is assumed to be in singly ionized form, as C^. We adopt 
the standard local value for the dust-to-gas ratio of 1:100 
(for further discussion see, e.g. Clover et al.||201Q ). We set 
the cosmic ray ioniz ation rate (CRIR) of atomic hydrogen 
to C = 3 X 10“^^ s“^ ( Clark et al.|2013 ). This value is a fac¬ 
tor ~ 1000 higher than the value measured in dense clouds 















































































3 


Model name 

Virial a 

Initial number density no 
[cm-3] 

3D velocity dispersion a„ 
[km/s] 

Radius R 

[pc] 

Free-fall time ffj 
[Myr] 

Simulation end tend 
[Myr] 

GC-0.5-1000 

0.5 

1000 

5.4 

8.9 

1.40 

1.00 

GC-2.0-1000 

2.0 

1000 

10.8 

8.9 

1.40 

1.40 

GC-8.0-1000 

8.0 

1000 

21.6 

8.9 

1.40 

1.40 


Table 1. Overview of the initial conditions for our different cloud models. 


in the solar neighbourhood (van der Tak & van Dishoeck 


2000) and is comparable to the high value that has previ¬ 


ously been inferred in the CMZ (Yusef-Zadeh et ah 2007 
Clark et al.||2013 ). Note, however, that it is only around a 


factor of 100 higher than the value derived in more diffuse 


gas in the solar neighbourhood (Indriolo et ah 2015). For 


the interstellar radiation field, we adopt the same spectral 


shape as given in Draine (1978). We denote the strength of 


the Draine ISRF as Go = 1 and perform simulations with 
a field strength Go = 1000 ( Clark et al.||2013 ). This corre¬ 
sponds to an integral flux in the energy range 6 — 13.6 eV of 


2.7 X 10 ^ erg cm ^s 


2.2 Model parameters 

The clouds that we model are initially spherical and located 
at the center of a large box of low-density gas. For the box, 
we use periodic boundary conditions. However, this is sim¬ 
ply a convenient choice, since the boundaries do not influ¬ 
ence the dense cloud evolution, because the size of the box is 
taken to be much larger than the size of the cloud. The cloud 
has a uniform initial hydrogen nuclei number density, which 
we set to 10^ cm“^, and a total mass of Mtot = 1.3 x 10^ M©. 
The initial cloud radius is R ^ 8.9 pc. In all of our sim¬ 
ulations, the density of the gas surrounding the cloud is 
1 cm“^, but we note that our results are insensitive to this 
value provided that it is much smaller than the mean cloud 
density. The cubic side length of the total simulation domain 
is set to 5x the individual cloud radius of 8.9 pc, i.e. 44.5 pc 
in total. All clouds have zero bulk velocity and are placed 
in the center of the box. The Voronoi cells initially have ap¬ 
proximately constant volumes. We initially start with 2x10® 
cells in total. The initial cell mass within the cloud corre¬ 
sponds to 2 Mq. 

We make use of a Jeans refinement criterion in order to 
accurately refine dense and collapsed gas regions in the box 
over the whole simulation time. We use a constant number of 
8 cells per Jeans length, which is sufficient to avoid artificial 


fragmentation (see, e.g. Truelove et al. 1998, Greif et al. 2011 
as well as Federrath et al. 2010b). Our code also includes a 
sink particle formation algorithm ( Bate et al.|1995 Jappsen 
et al. [20(55 Greif et al.|2011| ) to properly track the formation 
of stars during each run. Although our numerical models 
show active star formation by the end of the simulations, we 
ignore the sink particles in this study, since we only want to 
focus on the physical properties of the different gas phases. 
A more detailed analysis of the star formation history of 
our runs and the sink particle algorithm that we use can be 


found in Bertram et al. (2015). 


We assume that the gas has a turbulent velocity field 
with an initial power spectrum P{k) oc k~^ and consists 
of a natural mixture of solenoidal and compressive modes 


(see, e.g. Federrath et al.| 2010a). The strength of the tur¬ 
bulence is set by our choice of the initial virial a param¬ 
eter. This is defined as a = F^kin/|^pot|, i.e. the kinetic 
energy F^kin = l/2Mtot<7^ divided by the potential energy 
F^pot = —3GMt^ot/(5R) measured at the beginning of each 
run. The quantity Gv denotes the 3D rms velocity disper¬ 
sion in the cloud. We note that the virial parameter is also 
often calculated as a = 2F^kin/|F^pot| in the literature, which 
is different from the notation used in this paper (see, e.g. 
Federrath Klessen||201^ . With our definition, a value of 


a = 0.5 corresponds to a cloud in virial equilibrium and 
a > 1.0 describes clouds that are gravitationally unbound. 
In order to span a large range in the virial parameter space, 
we analyze models with a — 0.5, 2.0 and 8.0. Note that al¬ 
though it is convenient to parameterize the models in terms 
of a, the underlying physical quantity that changes as we 
change a is the turbulent velocity dispersion of the gas, since 
we keep Mtot and the initial radius R fixed in all three mod¬ 
els. It is therefore the changes to the velocity field, and the 
consequent changes to the density distribution of the gas, 
that dictate the observed trends seen in our simulations. 
However, for simplicity and clarity, we will simply refer to 
the different a models in our discussions in the following 
sections. 

The turbulence is not driven and hence decays through¬ 


out the simulation (Mac Low et al. 1998). The timesteps 
between the individual snapshots are Atgnap ~ 7 kyr, corre¬ 
sponding to ~ 200 snapshots per simulation in total. Table 
[T] summarizes the initial conditions of our clouds. We note 
that due to the high computational cost of our a = 0.5 
model, we had to stop our run already at t ~ 0.7 tff, while 
the other runs end at one free-fall time. This is due to an 
intense refinement of dense gas regions within the cloud. 
Furthermore, regarding the densities in the GMZ, we have 
to keep in mind that a significant number of GMZ clouds 
have densities higher than those modelled in our simula¬ 
tions, which would lead to even larger turbulent velocities 
in our simulations for the same constant total mass. 


2.3 Radiative transfer 

We post-process our data in order to generate synthetic 
line emission maps for the atomic and molecular tracers 
listed in Table using the radiative transfer code Radmc- 


3Ep| ( |Dullemond|2012 ). We use the Large Velocity Gradient 
(LVG) approximation (Sobolev 1957) to compute the G^, 
O and GO level populations. The LVG implementation in 
Radmc-3d is described in Shetty et al.| (2011a). We use the 
rate coefficients for collisional excitation and de-excitation 


1 www.ita.uni-heidelberg.de / ~dullemond / software / radmc-3d / 
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by atomic and molecular hydrogen tabulated in the Leiden 
database ( Schoier et al.||2005 ). 

At the moment, Radmc-3d cannot deal with Arepo 
data directly. Thus, we have to map the simulation out¬ 
put onto a cubic grid. For that, we adopt a resolution of 
512 cells in all three spatial dimensions, corresponding to a 
cell size of Ax = Ay = A; 2 ; = 0.087 pc. Moreover, we also 
use 512 channels in velocity space for our radiative transfer 
post-processing, corresponding to a spectral resolution of 
0.057 km s“^, 0.066 km s“^ and 0.102 km s“^ for our mod¬ 
els with a = 0.5, 2.0 and 8.0. Our choice of spatial resolution 
is determined by technical limitations within Radmc-3d: as 
it is a serial code, our grid size is limited by the requirement 
that the data should all fit within the memory of the com¬ 
puter that we use to carry our the post-processing. However, 
for most of our tracers, we find only minor differences in the 
emission as we increase the grid resolution. In particular, we 
find that a resolution of 512^ grid cells is enough in order to 
properly recover the emission of the cloud for all of the trac¬ 
ers apart from the [Ol] 63/im line (see Appendix]^. Hence, 
our main conclusions are not strongly affected by our choice 
of spatial resolution. 

Our simulations do not explicitly track the abundance 
of and so we need a procedure to relate the 

number density to that of ^^CO. A common assumption is 
that the ratio of to is identical to the elemen¬ 


tal abundance ratio of to (see, e.g. Roman-Duval 


et al. 2010). We make the same assumption and set the 
to ratio to a constant value, R 12/13 = 60. In 


reality, physical effects like chemical fractionation (Watson 
et al.|p'976 ) and photodissociation (Visser et al. 2009) will 


alter the abundances of relative to ^^CO, leading to 


a variable va lue of R 12/13 within the cloud. However, Sziics 
et al. (2014) have shown that the resulting maps do 


not differ greatly from those that we would obtain by as¬ 
suming a constant ratio of to ^^CO, and so we expect 

this approximation to be adequate. 

The radiative transfer calculation yields position- 
position-velocity (PPV) cubes of brightness temperatures 
Tb, which are related to the intensity via the Rayleigh-Jeans 
approximation, 

h 


r,M = (£)' 


2fcB’ 


( 2 ) 


where is the specific intensity at frequency v and ks the 
Boltzmann constant. For each run, we carry out the post¬ 
processing for each of the tracers listed in Table along an 
arbitrary line of sight (LoS), using the last snapshot of each 
model. Since we use isotropic turbulence without magnetic 
fields, we do not expect our results to significantly depend 
on the specific choice of the LoS and hence we only focus on 
emission observed along the ^-direction (see also Appendix 

H- 


3 RESULTS 

In this section, we analyze various physical cloud parame¬ 
ters using the latest time snapshot (see Table . We start 
with the thermal state of the MCs (Section [3.1[ ) and investi¬ 
gate which chemical components best trace the total shape 
of the cloud (Section 3.2). Afterwards, we estimate the ef¬ 
fective cloud radii by using synthetic observations in various 


Tracer 

Type 

Transition 

A 

[fim] 

V 

[GHz] 

12CO 

Molecular 

J =1^0 

2600 

115 

13CO 

Molecular 

J =1^0 

2720 

110 

[CII] 

Atomic 

^P3/2 —1/2 

158 

1900 

[Oi] 

Atomic 

3Pi 

63 

4744 

[Oi] 

Atomic 

3Po 

145 

2060 


Table 2. Tracer name and type, quantum mechanical transition, 
wavelength A and frequency u of various fine structure lines, which 
we model with the radiative transfer code Radmc-3d. 


diagnostic lines (Section |3.3| ). We then explore the ability of 
the lines listed in Table | 2 | to trace dense MC regions (Sec¬ 
tion |3.4| and continue with a quantitative estimate of the 
mass fraction of the cloud that is traced by the line emission 
(Section 3.5). At the end, we discuss the Aco-factor of the 
MCs. 


3.1 Thermal state of the clouds 


In Fig. we show the two-dimensional PDFs of gas tem¬ 
peratures and densities within the clouds for the different 
virial parameters and for all chemical components: total 
mass, molecular hydrogen, atomic oxygen, ionized carbon 
and carbon monoxide. These PDFs are mass-weighted, with 
the colour-coding indicating the fraction of the total mass 
of the cloud (in the upper row of panels) or of the species 
of interest (in the remaining panels) located at each point 
in density-temperature space. The number densities shown 
in the upper row of panels are the total number density of 
particles; in the remaining panels, we show instead the num¬ 
ber density of the species of interest no, etc.). Figure 

[ 1 ] demonstrates that the thermal state of the clouds for all 
of the components only weakly depends on the virial pa¬ 
rameter. However, we find strong differences in the PDFs 
between the various tracers. 

We find a significant amount of warm gas at densities 
n < 1000 cm“^ and temperatures between T ^ 300 — 10^ K. 
We also note that there is a substantial scatter in the tem¬ 
perature at every value of the density in this regime. This 
gas corresponds to the warm and tenuous envelope in the 
outer regions of the clouds with low extinction. It is heated 
by the strong external ISRF via the photoelectric effect. For 
densities higher than n > 1000 cm“^, we find an inverse 
relationship between the temperature and the density. The 
temperature decreases with increasing density, owing to the 
efficient self-shielding of the dense gas from the external ra¬ 
diation field. Such high densities only occur in the inner re¬ 
gions of the clouds with temperatures of T < 100 K, which 
are cold enough to result in active star formation (for further 
information about the star formation history of our numeri¬ 


cal simulations, we refer the reader to Bertram et al. 2015). 
These are also the parts of the cloud in which we find molec¬ 
ular H 2 and CO gas in our simulations. 

If we compare the thermal state of the total mass to 
the thermal state of the other tracers, we find that only the 
atomic oxygen and ionized carbon show a similar behavior 
over the whole range of temperatures. This is because these 
components exist in the diffuse cloud regions as well as in 
the denser, more shielded gas. We also see that molecular 
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Figure 3. Total mass per unit logarithmic column density for 
the total and the H 2 gas as shown in Fig. I^for our different virial 
parameters: a = 0.5 (top), 2.0 (middle) and 8.0 (bottom). 


hydrogen is strongly dissociated in the diffuse, warm enve¬ 
lope of the cloud. It only reaches maximal temperatures of 
~ 4000 K, which is lower than the warmest regions observed 
for the total gas. In the case of CO, the range of temper¬ 
atures (30 — 300 K) is even smaller. This suggests that CO 
traces a denser, colder part of the cloud than [Cii] or [Oi]. 
We will examine the consequences of this in more detail in 
the following sections. 


that the H 2 column density traces the total column density 
fairly well. In the runs with a = 2 and a = 8, however, 
the cloud is not gravitationally bound and hence starts to 
expand. In addition, the strong turbulence opens up large 
channels into the cloud by lowering the column density lo¬ 
cally, allowing radiation to penetrate deeper into the MC, 
leading to enhanced photodissociation of the molecular hy¬ 
drogen. Hence, we do not find a good correlation between the 
total and the H 2 column density in the a = 2.0 and a = 8.0 
runs. Molecular hydrogen can only exist in the most dense 
regions of the cloud, where it is best able to self-shield from 
the strong external radiation field. 

This is also illustrated in Fig. which shows the cor¬ 
responding column density PDFs for three virial parame¬ 
ter models. If we compare the H 2 column density PDFs 
with each other, we find that the peak column density in 
the a = 0.5 model is shifted to larger values compared 
to the a = 8.0 model by about one order of magnitude. 
This is because the turbulent kinetic energy in the virialized 
Q/ = 0.5 run is much smaller than in the a = 8.0 run, lead¬ 
ing to denser cloud regions and thus to larger average col¬ 
umn densities. Moreover, the fraction of H 2 gas in the total 
box decreases with increasing virial parameter, owing to the 
high external ISRF, which photodissociates those molecular 
components more effectively as the radiation can penetrate 
deeper into the cloud’s interior. Furthermore, we see that 
the H 2 gas generally extends to much lower column densi¬ 
ties than the total gas. This is because the total gas consists 
of both atomic and molecular hydrogen and thus does not 
fall below lower values than A^tot ~ 10^^ cm“^ due to the Hi 
background, while the fraction of H 2 gas can continuously 
decrease in the outer regions due to photodissociation. 

In analogy to Fig. Fig. shows velocity-integrated 
intensity maps computed along the LoS in the z-direction 
for the same virial parameters and the different tracers pre¬ 
sented in Table We again find that the molecular tracers 
and are photodissociated at the edges of the MC 

by the strong external ISRF. This effect is even stronger if 
the amount of turbulent kinetic energy in the box is high. 
Comparing the different tracers to the total gas column in 
Fig. we find that the cloud’s shape is well reproduced by 
all atomic tracers, although the integrated brightness tem¬ 
peratures of these lines are significantly smaller than those 
of the ^^CO and ^^CO lines. This holds for all three virial 
parameter models. Conversely, in the case of CO, we find 
that neither of its isotopologues describes the total gas dis¬ 
tribution of the cloud. However, we observe that they do 
remain reasonably good tracers of the dense H 2 , owing to 
the effect of self-shielding. We will quantify this finding in 
Section 13.51 


3.2 Tracing the shape of the MC 

Fig. shows logarithmic column density maps for the to¬ 
tal density and the H 2 density computed along the LoS in 
the ^-direction for virial parameters a = 0.5, 2.0 and 8.0 for 
the last time snapshots (see Table . We see that as we 
increase o, and hence the velocity dispersion of the gas, the 
structure of the cloud changes significantly. When o = 0.5 
and the cloud is gravitationally bound, it remains relatively 
compact. It is therefore able to shield itself relatively well 
from the effects of the high external ISRF, with the result 


3.3 Estimating the effective cloud radii 

In the previous section we established that carbon monoxide 
is not a good tracer for the total column density, while the 
atomic species better reflect the distribution of the total gas 
mass. In this context, CO tends to signihcantly underesti¬ 
mate the radius of the total cloud in such an extreme envi¬ 
ronment. This can compromise any observational estimates 
of the virial parameter a, which requires an accurate esti¬ 
mate for the cloud’s radius (see also our follow-up study in 
Paper H). However, we note that it is generally complicated 
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Figure 1. Physical state of the clouds at the end of our simulations. We show two-dimensional PDFs of gas temperature and number 
density for various chemical components for runs with a = 0.5 (left column), 2.0 (middle column), and 8.0 (right column). In the upper 
row, the number densities shown are the total particle number density. In the remaining rows, on the other hand, they are the number 
density of the chemical species indicated in the panel (e.g. etc.). Note that the scaling of the density axis therefore differs from 

panel to panel. The color coding shows the fraction of the total gas mass in each logarithmic density and temperature bin. 
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a = 0.5 a = 2.0 



a = 8.0 



Figure 2. Logarithmic column density maps computed along the LoS in the ^-direction for the total density (top row) and the H 2 
density (bottom row) for different virial parameters: a = 0.5 (left column), 2.0 (middle column) and 8.0 (right column). Each side has a 
length of 44.5 pc. The initial cloud radius of 8.9 pc is indicated in the bottom panels. 


Model 

[Cll] (158/^m) 
[Kkms-i] 

[Ol] (145/im) 
[Kkms-i] 

[Ol] (63/im) 
[Kkms-i] 

(2600 ^im) 

[Kkms-i] 

(2720/im) 
[Kkms-i] 

GC-0.5-1000 

3.0 

1.0 

0.6 

1.0 

1.0 

GC-2.0-1000 

3.0 

1.0 

0.6 

0.3 

0.3 

GC-8.0-1000 

3.0 

0.8 

0.5 

0.1 

0.1 


Table 3. Threshold values Wthresh used for the integrated intensity maps in Fig. ^ in order to estimate the effective cloud radius 
The values are chosen such that the denser cloud regions can be clearly separated from the diffuse gas phase around it, leading to a 
reasonable estimate for the effective cloud radii presented in Table ^ 


to define a proper cloud radius. Molecular clouds are com¬ 
plex hierarchical systems and estimates of cloud radii always 
come with some assumptions of how to define outer cloud 
boundaries. Hence, if we refer to the effective cloud radius 
in the following, we always mean the radius of a spherical 
cloud with the same surface area. 

Below, we try to estimate the effective cloud radius i^eff 
in a similar fashion to what would be applied to observa¬ 
tional data sets. At first sight, this seems to be complicated, 
since all clouds shown in Fig.^are far away from a spherical 
geometry. Nevertheless, we can estimate an effective radius 
by adopting an integrated intensity threshold Wthresh for 
each map and counting the number of cells N, which fulfill 
Wcioud > Wthresh- The valucs for Wthresh are listed in Table 
l^and are chosen so that there is a clear separation between 
emission from the cloud and from the diffuse gas in which it 
is embedded. Our results are not particularly sensitive to the 
exact value chosen for Wthresh, as long as it is large enough 
to not be contaminated by the diffuse emission, but not so 
large that we miss significant emission from the cloud itself. 


For example, in the case of model GC-0.5-1000, values of 
Wthresh anywhere in the range 2 < Wthresh < lOKkrns”^ 
give essentially the same results for For comparison, 

we also give the effective radii of the clouds that incorpo¬ 
rate ~ 80% of the total gas mass. 

We have also verified that our adapted thresholds are 
comparable or greater than the practical thresholds that we 
would obtain using modern instruments to measure the line 
emission. For example, in a single pointing, the GREAT in¬ 
strument (German REceiver for Astronomy at Terahertz 
Frequencies, Heyminck et al. 2012) on-board the Strato¬ 
spheric Observatory for Infrared Astronomy (SOFIA) can 
measure integrated intensities of the [Gll] line down to 
around ~ 0.06Kkms“^ in around one hour of integration 
tim^ although multipixel array detectors are available to¬ 
day that become more and more efficient. Mapping an entire 
cloud using multiple pointings obviously allows one to spend 
much less time per pointing, yielding a higher threshold, but 


^ https://great.sofia.usra.edu/cgi-bin/great/great.cgi 


([ LUO] / ([ LUO] / 2Hn)6o| 
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Model 

80% of total mass 

[pc] 

[Gii] (158 /xm) 

[pc] 

[Ol] (145/im) 

[pc] 

[Ol] (63 /xm) 

[pc] 

(2600(Um) 

[pc] 

i3CO (2720/im) 

[pc] 

GC-0.5-1000 

11.0 

9.0 

8.7 

8.1 

6.6 

6.2 

GC-2.0-1000 

13.1 

13.8 

13.0 

12.2 

7.6 

6.8 

GG-8.0-1000 

17.5 

20.1 

19.9 

19.7 

8.4 

6.9 


Table 4. Effective radii of the clouds shown in Fig.^ as inferred from equation]^ based on the threshold intensities given in Table 
For comparison, we also give the effective radii of the clouds that incorporate ~ 80% of the total gas mass in the second column. 


it remains plausible to map extended regions with a sensi¬ 
tivity comparable to or better than the adopted threshold. 

We compute the effective cloud radius i^eff by assuming 
an effective cloud area Agff = where Agff = NAx^ and 

Ax is the cell size introduced in Section [2.3| Thus, we get 



Table [4] shows the different effective radii of the various 
clouds presented in Fig. computed via equation (§. If 
we compare the different values listed in Table we find 
that the cloud radii measured with the integrated intensi¬ 
ties of and are significantly smaller by a factor 

of 1.5 — 2.5 compared to the values measured with the in¬ 
tegrated intensities of the atomic tracers. This is because 
the strong ISRF and high CRF destroy the molecular trac¬ 
ers in the more diffuse regions, leading to smaller effective 
radii. Conversely, the radii measured with the atomic com¬ 
ponents roughly agree with each other for one specific virial 
parameter model and give the best radius estimate for the 
total cloud that incorporates ~ 80% of all gas. These values 
confirm our assumption that the atomic species accurately 
reflect the spatial distribution of the total gas column and 
that they yield a robust estimate of the effective radius of 
the cloud. 


3.4 Tracing the total and H 2 column density 

Now we turn to the densest parts of the cloud, because this 
is where stars are formed, and we study how well atomic 
fine-structure lines can be used to trace this regime. Fig. 
andj^show velocity-integrated intensities plotted against the 
total column density and H 2 column density for the different 
virial parameters and all tracers presented in Table In 
these plots, we also show the the mean value at each column 
density, which is indicated by a black solid line. The small 
fluctuations at high column densities in each plot are caused 
by the small number of cells available to compute an average 
value of W. 

We find clear differences for the various chemical com¬ 
ponents. For the total density in Fig. we observe that 
the and tracers fall off sharply with decreasing 

column density below ~ lO^^cm”^, which is due to pho- 
todissociation of the gas in the diffuse cloud regions, owing 
to the strong external radiation field. On the other hand, 
at high column densities, we find a saturation of the 
emission due to the effect of the line opacity, while the 
emission is a significantly better tracer of the column density 
than ^^CO. This is because is optically thin and con¬ 

sequently can better trace compact regions within the cloud. 
However, even in a small regime where Atot and Wco show 


some degree of correlation, we find a substantial amount of 
scatter around the mean value for both and ^^CO. 

Furthermore, we also show the scatter plots for the different 
atomic species in Fig.[^ As expected, we see that all atomic 
components extend to significantly lower column densities 
of ~ 10^^ cm“^ compared to carbon monoxide. Hence, both 
atomic carbon and oxygen are much better tracers of the 
low column density material of the total gas in the cloud. 
However, both the [Cll] 158 /xm and the [Of] 63 /xm line sat¬ 
urate at higher column densities, owing to the effect of the 
line opacity and thus either or the [Of] 145/xm line is 

to be preferred as tracers for the dense cloud material. 

For the H 2 column density in Fig. we find a differ¬ 
ent result. Firstly, we see that both and extend 

to significantly lower H 2 column densities of ~ 5 x 10^° — 
10^^ cm“^, reflecting the fact that CO is a better tracer of 
molecular mass than of total mass. At higher column densi¬ 
ties, we find a qualitatively similar behavior as for the total 
gas in Fig. again owing to the effect of the line opacities. 
Moreover, we also show the scatter plots for the different 
atomic species. We note that both the atomic oxygen and 
carbon extend over a much larger range in H 2 column den¬ 
sity than in total gas column density, showing that those 
components can also be used in order to accurately trace 
a large fraction of the molecular hydrogen mass. However, 
the correlation between the emission of [Cll] 158 /xm, [Of] 
63 /xm and Ahs breaks down at H 2 column densities above 
~ 10^° — 10^^ cm“^, owing to the optical depth of these lines 
in this density regime. This means that one has to rely on 
other probes in order to infer information about the com¬ 
pact gas regions. On the contrary, the correlation between 
the [Of] 145 /xm emission and seems to be provide a 
better way to study the dense H 2 clumps in the cloud. 

We also analyze the impact that our over-estimation of 
the CO abundance may have on the relationship between the 
integrated intensity of the atomic tracers and the column 
density (see Appendix . However, the influence of this 
turns out to be small. 


3.5 Mass fraction traced by the emission 

In Section |3.2| we speculated that a significant fraction of 
the total gas mass may not be traced by ^^CO or ^^CO. 
To quantify this finding. Fig. shows the fraction of total 
and H 2 gas mass observed along the LoS having velocity- 
integrated intensities greater than a minimum value Wmin 
for our various tracers. Furthermore, we analyze how this 
mass fraction changes as we vary Wmin- In this context. Fig. 
[^defines a theoretical limit for detecting gas above a given 
sensitivity threshold, e.g. determined by the telescope or the 
detector. 
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Figure 4. Velocity-integrated intensity maps computed along the LoS in the 2 ;-direction for different virial parameters: a = 0.5 (left 
column), 2.0 (middle column) and 8.0 (right column). From top to bottom: integrated intensity maps for (J = 1 —)■ 0), 

(J = 1 —)■ 0), [Ol] (63/rm), [Ol] (145/rm) and [Cll] (158p,m). Note the different scaling in the colorbars at the right hand side. Each side 
has a length of 44.5 pc. The initial cloud radius of 8.9 pc is indicated in the bottom panels. The contour lines show the threshold values 
given in Table 
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Figure 5. Log-log plot showing velocity-integrated intensities against total column density for different virial parameters: a = 0.5 (left 
column), 2.0 (middle column) and 8.0 (right column). Shown are all our tracers presented in TableIn each plot, the mean value is 
indicated by a black solid line. 
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Figure 6. Same as Fig. but with the column density of H 2 . 
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Figure 7. Cumulative fraction of the total (top row) and H 2 (bottom row) mass in the cloud for different virial parameters a = 0.5 (left 
column), 2.0 (middle column) and 8.0 (right column), traced by our various chemical components presented in Tableas a function of 
the minimum velocity-integrated intensity in the line. 


We see that all lines are a strong function of the min¬ 
imum velocity-integrated intensity and of the virial param¬ 
eter. In case of the mass fraction of the total gas, we find 
that our different atomic tracers recover almost all of the 
total mass if the minimum integrated intensity is Wmin ~ 
0.1 Kkms“^ in the a = 0.5 model and Wmin ~ 1.0Kkms“^ 
in the a = 2.0 and a = 8.0 models. Furthermore, as already 
indicated in Section [3^ we see that in the case of very small 
threshold values Wmin in the a = 0.5 model, both and 
only trace about 75% of the total mass or less. This 
fraction decreases even further for higher a values. For ex¬ 
ample, if we adopt Wmin ~ O.lKkms”^ in the cloud with 
a = 8.0, we find that only ~ 40% of the total mass is traced 
by CO, while all atomic tracers recover almost 100% of the 
total gas mass. Furthermore, we observe narrow emission 
ranges for some of the atomic components, which trace al¬ 
most 100% of the total mass. For example, regarding our 
intermediate a = 2.0 model and the [Ol] 63 /xm tracer, al¬ 
most 100% of the total mass lies in the narrow emission 
range between 0.5 < Wmin < 3.0. In the case of the other a 
models, we find a similar behavior. This is because there are 
almost no lines-of-sight having higher values of integrated 
intensity (see also Fig. and [^, which may be different in 
even denser clouds. 

Regarding the H 2 mass, we see that both and 

trace a higher H 2 mass fraction compared to the total 
mass of the cloud. This is because the molecular hydrogen 
can only exist in dense regions of the cloud, owing to the 
ability to self-shield from the external radiation field. This is 
similar to our finding for the carbon monoxide. Furthermore, 
we observe that the atomic components recover almost 100% 


of the H 2 mass if we set Wmin ~ 0.5Kkms“^. Hence, we 
find that those components again trace a significantly higher 
mass fraction of the gas compared to the carbon monoxide. 
Regarding both and ^^CO, we see that the amount 

of mass recovered by those molecular tracers decreases with 
higher virial parameter for a fixed minimum value of the 
integrated intensity. For example, if we adopt a value of 
Wmin ~ O.lKkms”^, almost 98% of the H 2 mass is traced 
by for our a = 0.5 model, while ~ 90% is traced in 

our a = 2.0 and only ~ 70% can be recovered in the a = 8.0 
model. Moreover, similar to the mass fraction of the total gas 
mass, we also find a narrow emission range for the various 
atomic lines in all a models, which trace about 100% of the 
H 2 mass. In contrast, the emission of the and 

molecules is extended over a much wider range in integrated 
intensities. Hence, a much larger emission range is needed 
in order to trace the whole mass of molecular hydrogen in 
the cloud. 

In summary, we conclude that CO is only a good tracer 
for the molecular content of the cloud, but can miss a signifi¬ 
cant fraction of the total gas mass, particularly for unbound 
clouds. 


3.6 Estimating the CO-to-H 2 conversion factor 

Based on our maps of column densities for H 2 (Fig. and 
velocity-integrated intensities for (Fig. |^, we also es¬ 

timate the value of the CO-to-H 2 conversion factor, Xco, 
in our various models. Xco is a widely used quantity in 
particular in extragalactic astronomy to derive H 2 column 
densities from CO observations (Bolatto et al. 2013). It is 
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Figure 8. Logarithmic maps of the CO-to-H 2 conversion factor, 
= ^H 2 /^C 0 : for fho three models: a = 0.5 (top), 2.0 
(middle) and 8.0 (bottom). The average XcQ-factor values are 
given in Table We see that only the inner parts of the cloud 
yield the canonical Xco“fcictor value on average, while the outer 
cloud regions exhibit significantly larger Xco values. The gray 
background denotes regions where the CO emission is zero. Hence, 
it is impossible to compute a value for Xco there. 




defined via 




VLco,i-o 


(4) 


where iVH 2 denotes the H 2 column density and Wco,i-o the 
velocity-integrated brightness-temperature. Previous stud¬ 
ies tried to estimate the Xco-factor in both observations 


and numerical simulations (see, e.g. [Solomon et al. 

1987 

Polk et aL|1988| [Young X Scoville[[1991[ Dame et al. 

2001 

Liszt et aL[2010| [Leroy et aL|2011[ [Glover X Mac Low 

2011 


Model 

Xco 

[cm-^K-lkm-ls] 

GC-0.5-1000 

3.9 X 1020 

GC-2.0-1000 

2.6 X lO^o 

GC-8.0-1000 

1.3 X 1020 


Table 5. Values of the CO-to-H 2 conversion factor, relating the 
amount of (J = 1 ^ 0) emission Wqq to the H 2 col¬ 

umn density for all our a runs. For the computation of the 
Vco“factor, we only include pixels where the velocity-integrated 
intensity for carbon monoxide is larger than the corresponding 
threshold given in Table Thus, Xco only accounts for regions 
with a significant amount of CO emission in the cloud. 


Glover Clark||2012| [Bolatto et al.|[20T3l [Narayanan <V 


Hopkins||2013| [Clark L Glover||2015f ^^ For example, [Shetty 


et al.[ (2011a,b) investigated the Xco-factor in numerical 


simulations of turbulent clouds, varying environmental prop¬ 
erties such as the initial number density, the metallicity 
and the external UV field, finding cloud average values 


Xco ~ 2 — 4 X 10^^ cm“^ km“^ s for solar-metallicity 
models. Moreover, Bolatto et al. (2013) review the eflPorts 
to measure the CO-to-H 2 conversion factor for molecular 
clouds in the Milky Way disk, recommending a similar factor 
Xco ~ 2 X 10^° cm“^ km“^ s with ±30% uncertainty. 
However, all these studies arrived at the conclusion that the 
Xco-factor in the Galaxy is remarkably constant, having a 
value of Xco ~ 2 — 4 x 10^° cm“^ km“^ s. 

Nevertheless, in the last few years concerns have been 
raised that the Xco-factor may be significantly different 
from the canonical value in MCs located in the galactic 
center. For example, Narayanan et al. (2012) have shown 
that large gas temperatures and velocity dispersions can in¬ 
crease the GO intensity, thus decreasing Xco- On the other 
hand, the strong ISRF and GRF can also destroy the car¬ 
bon monoxide, increasing Xco (Glark X Glover 2015). In 
addition, observational measurements of Xco at the cen¬ 
ter of other nearby spiral galaxies typically find values that 


are lower than the canonical Galactic value (Sandstrom 


et al. 2013). Moreover, conditions in the Galactic Genter 


are somewhat similar to those found in ULIRGs, and there 
is considerable observational evidence that Xco is smaller 
in ULIRGs than in normal spiral galaxies (see e.g. the de¬ 
tailed discussion in Bolatto et al.[[2013 ). So far, we are still 
missing a reliable picture of the Xco-factor in different (and 
extreme) physical environments. 

In the following, we evaluate the Xco-factor for our 
three a models and list our results in Table In general, 
the Xco-factors given in Tableare computed by taking the 
ratio of the mean values of the H 2 column density and the 
velocity-integrated intensity of GO. We only include pixels 
where the velocity-integrated intensity for carbon monoxide 
is larger than the corresponding threshold given in Table 
These are 1.0, 0.3 and 0.1Kkms“^ for the models with 
a = 0.5, 2.0 and 8.0. Thus, Xco only accounts for regions 
with a significant available amount of GO emission in the 
cloud, as illustrated in Fig. 

We find a range of values for the GO-to-H 2 con¬ 
version factor in Table reaching from ~ 1 — 4 x 

10^° cm“^ km“^ s. Our estimates are in good agreement 
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to estimates of the canonical Xco-factor in the Milky Way. 
However, we again emphasize that our [Nelson Lange^ 
( |l997l l network overestimates the rate at which CO forms 
(see also Section and hence we expect that our Xco- 
factors listed in Tableware lower limits. 

Furthermore, we show logarithmic maps of the Xco- 
factor for all our a models in Fig. Here, the Xco-values 
are evaluated by computing the ratio of the H 2 column den¬ 
sity and the velocity-integrated intensity of CO for each in¬ 
dividual pixel. The gray background denotes regions where 
the CO emission is zero. Hence, it is impossible to compute 
a value for Xqo there. We find significantly different values 
for the inner and the outer parts of the clouds. While the in¬ 
ner cloud regions reproduce the canonical Xco-factor value 
on average, we find that the outer cloud regions exhibit sig¬ 
nificantly larger Xco values compared to the canonical one. 
This is because the carbon monoxide is photodissociated at 
the edges of our clouds by the strong ISRF, leading to an 
increase in Xco- In contrast, the molecular CO is better 
able to self-shield in the dense parts of the cloud and hence 
we observe that those regions fairly reproduce the canonical 
value. 


4 SUMMARY AND CONCLUSIONS 


In this paper we have analyzed synthetic images of MCs in a 
CMZ-like environment. For this purpose, we have performed 
numerical simulations of model clouds with the moving mesh 
code Arepo ( Springel|2010 ) using environmental properties 
comparable to those experienced by typical CMZ clouds. We 
adopted values for the interstellar radiation field (ISRF) and 
the cosmic-ray flux (CRF) that are a factor of ~ 1000 larger 
than the values measured in the solar neighbourhood ( [Clark 
et al. 2013). We simulated clouds with an initial number den¬ 
sity of no = 10^ cm“^ and studied the impact of different 
virial a parameters of a = 0.5, 2.0 and 8.0. The total cloud 
mass was set to a constant value of Mtot = 1.3 x 10^ M©. 
Furthermore, we used the radiative transfer code Radmc-3d 
( Dullemond|2012| to compute synthetic maps of our clouds 
in important diagnostic lines, i.e. we model the emission 
of the cloud in [Cll] (158/im), [Ol] (145/xm), [Ol] (63/xm), 
^^CO (2600 /xm) and ^^CO (2720 /xm). We report the follow¬ 
ing findings: 

• Atomic carbon and oxygen is found in both cold and 
warm regions in the clouds, but CO only traces cold regions 


(see Section 3.1). 

• If the cloud is virialized {a = 0.5), the H 2 gas is much 
better able to self-shield from the high external ISRF com¬ 
pared to other runs with a larger a value. This is because 
the gravitational pressure compresses the gas to much higher 
densities because the amount of turbulent kinetic energy in 
the box is smaller in the low a models than in the high 
virial parameter models. In the latter case, radiation can 
penetrate deeper into the MC, leading to photodissociation 


of the molecular hydrogen (see Section 3.2). 


• We show that the atomic components trace the shape 
of the total cloud very well even in diffuse cloud regions (see 
Section [3.2[) . 


• On the other hand, the molecular gas (H 2 as well as CO) 
is photodissociated at the edges of the MC by the strong 
external ISRF. 


• As a consequence, the photodissociation has a strong 
impact on measurements of the effective MC radius using 
a molecular tracer, which tend to significantly underesti¬ 
mate the radius of the total cloud in environments where 
the external radiation field and the turbulent kinetic energy 
is high. However, we find that the atomic tracers recover 
a significantly larger fraction of the total cloud mass (see 
Section 3.3). 

• We find a large dynamical range for the [Oi] 145 /xm 
line in column density space, which traces both the H 2 and 
the total gas up to the most dense regions in the cloud. 
In contrast, the [Oi] 63/xm and [Cll] 158/xm lines already 
saturate at signihcantly lower column densities. Hence, we 
would not expect that those components trace the compact 
MC regions very well (see Section 3.4). 

• We find that all atomic components trace almost 100% 
of the H 2 and the total gas mass above a sensitivity threshold 


of Wr. 


O.lKkms ^ for the velocity-integrated intensi¬ 


ties. However, we also find that CO only traces a significantly 
lower mass fraction of H 2 and the total density compared to 
the atomic components (see Section 3.5). 

• We compute values of the CO-to-H 2 conversion factor 
= Nu^/Wco, which relates the H 2 column density iVH 2 
to the amount of ^^CO (J = 1 ^ 0) emission Wco- We find 
values in the range Aco ~ 1 — 4 x 10^° cm“^ km“^ s, 
in good agreement to the canonical Aco-factor values ob¬ 
tained in observations of MCs in the Milky Way (see Section 


3.6). However, we caution that simplifications in our chem¬ 


ical model mean that these values are lower limits. 
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APPENDIX A: RADIATIVE TRANSFER 
POST-PROCESSING: INFLUENCE OF GRID 
RESOLUTION 


Since Radmc-3d cannot deal with Arepo data directly, we 
have to map the simulation output onto a cubic grid. The 
results presented in this paper are all based on a grid reso¬ 
lution of 512^ cells. However, we examine the sensitivity of 
our results to the choice of the number of grid cells in each 
dimension. Therefore, we run a similar radiative transfer 
post-processing with Radm c-3d , using resolutions of 256^ 
and 128^ grid zones. In Fig. 


Al 


we compare the probability 
density functions (PDF) of the velocity-integrated intensi¬ 
ties for our various tracers for the different grid resolutions, 
while using the intermediate a = 2.0 cloud models as an ex¬ 
ample. We hnd that a resolution of 512^ grid cells is enough 
in order to properly recover the emission of the cloud for all 
of the tracers apart from the [Ol] 63/im line. In the case of 
[Ol] 63/im, the maps are well-converged for integrated inten¬ 
sities of below 1.5 K km s“^, but at higher integrated intensi¬ 
ties, increasing the resolution depresses the PDF. However, 
as we have already seen in Section [3.4| integrated intensities 
of this magnitude are only recovered for [Ol] 63/im along 
lines of sight where the line is already optically thick. The 
precise values of the integrated intensity that we recover in 
this regime therefore do not significantly affect our conclu¬ 
sions. 

For all of the tracers, we find that the differences be¬ 
tween the PDFs we recover from the 256^ and 512^ runs 
are always significantly smaller than the differences between 
the 128^ and 256^ models. Thus, we do not expect to find 
larger differences in the PDFs for models with even higher 
resolution, e.g. between 512^ and 1024^. Hence, due to the 
significant increase of the computational costs for the radia¬ 
tive transfer post-processing for even higher resolution runs, 
we focus on models with a resolution of 512^ grid cells in 
this paper, which is enough to demonstrate our basic con¬ 
clusions. 


APPENDIX B: DEPENDENCE ON VIEWING 
ANGLE 

When constructing our synthetic emission maps, we choose 
to focus on a line of sight (LoS) parallel to the z-axis of the 
simulation volume. However, as the turbulence in our sim¬ 
ulations is isotropic, we expect our results to be insensitive 
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Figure Al. Emission PDFs for our various tracers for different 
resolution models with 128^, 256^ and 512^ grid cells for clouds 
with a = 2.0. 
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Figure Bl. Example of emission PDFs for the tracer for 

the resolution model with 512^ grid cells. We choose the inter¬ 
mediate cloud model with ol — 2.0. The above figure shows the 
emission PDF for the x-^ y- and ^-direction. Although there are 
slight variations in the different PDFs, the bulk emission is the 
same for all directions. 


to this choice. To test this, we have also made maps of 
for LoS parallel to the x and y axes, using our standard grid 
resolution of 512^ zones, for the run with a = 2.0. The re¬ 
sulting integrated intensity PDFs are shown in Figure [BT] 
Although there are slight variations in the different PDFs 
for y and z, the bulk emission is the same for all direc¬ 
tions. We find a similar behavior also for the other atomic 
tracers. Hence, we conclude that it is enough to focus on one 
specific LoS, for example the ^-direction, in order to get an 
idea about the underlying physical parameters of the cloud. 
We note, however, that for magnetized clouds this is not 
necessarily true, as in this case the turbulence will no longer 
be isotropic, unless the field is very weak. 


APPENDIX C: INFLUENCE OF THE CO 
ABUNDANCE ON THE LINE EMISSION OF 
[CII] AND [OI] 


In Section m we emphasize that the simplified [Nelson <S] 
Danger ( |1997 ) network that we use in our study is known 


to somewhat overestimate the rate at which CO forms in 
our simulations. Thus, we may expect that this has a conse¬ 
quence on the distribution of the line emission of [Cll] and 
[Ol]. In order to quantify this, we perform additional runs of 
radiative transfer post-processing with updated abundances 
of carbon and oxygen for comparison by assuming zero CO 
abundances in our simulations. As an example, we plot the 
different line emission PDFs for our various atomic trac¬ 
ers for a fixed resolution of 512^ grid cells for clouds with 
a — 2.0. The result is shown in Fig. jClj We find that the 
[Cll] 158 /im and [Ol] 63 /im line emission are unaffected by 
the changed abundances. This suggests that in the region 
where the CO abundance is uncertain, these two lines are 
already optically thick and hence are in sensitive to this 
uncertainty. Support for this interpretation comes from the 
behavior of the [Ol] 145 /im line. As we have already seen, 
this has a much smaller optical depth, and hence is more 
strongly altered by the changed abundances. However, even 
in this case, the effect is relatively small: the brightest re- 





Figure Cl. Emission PDFs for our various atomic tracers for 
a fixed resolution of 512^ grid cells for clouds with a = 2.0. 
Shown are the PDFs given in Fig. |Al| (with CO) together with 
PDFs computed from additional runs of radiative transfer post¬ 
processing with updated abundances of carbon and oxygen as¬ 
suming zero CO abundances in our simulations (without CO). 
We do not find significant changes in the different emission PDFs 
of the atomic tracers. 


gions become around 0.1 dex brighter, but the majority of 
the PDF remains unaffected. We can therefore conclude that 
the inaccuracy introduced into our predicted [Cll] and [Ol] 
maps by our simplified chemical model is unimportant. 

















